Normalized entropy density of the 3D 3-state Potts model 
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Using a multicanonical Metropolis algorithm we have performed Monte Carlo simulations of the 
3D 3-state Potts model on L 3 lattices with L — 20, 30, 40, 50. Covering a range of inverse tem- 
peratures from /3 m in = to /3 m ax = 0.33 we calculated the infinite volume limit of the entropy 
density s(/3) with its normalization obtained from s(0) = In 3. At the transition temperature the 
entropy and energy endpoints in the ordered and disordered phase are estimated employing a novel 
reweighting procedure. We also evaluate the transition temperature and the order-disorder interface 
tension. The latter estimate increases when capillary waves are taken into account. 
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PACS numbers: PACS: 05.50.+q, 11.15.Ha, 12.38.Gc, 25.75.-q, 25.75.Nq 



I. INTRODUCTION 

The 3D 3-state Potts model plays a role in our under- 
standing of the properties of the deconfming phase tran- 
sition and the structure of QCD. The high-temperature 
vacuum of QCD is characterized by ordered Polyakov 
loops which behave similarly as spins in the low temper- 
ature phase of the 3D 3-state Potts model. The Polyakov 
loop serves as an order parameter of SU(3) pure gauge 
theory with its symmetry determined by the Z% center 
of the gauge group, and the phase transition is weakly 
first order. This maps naturally onto the order param- 
eter (magnetization) of the 3D 3-state Potts model 
For this and other reasons a number of numerical studies 
of the model were performed @, & II H & . 

In this paper we use the multicanonical technique [8| 
to calculate the properly normalized entropy all the way 
from (3 = to across the phase transition. Besides this 
new result a number of additional estimates (see the ab- 
stract), arc obtained and compared with the literature. 
Our paper is organized as follows: In section [TT] we intro- 
duce basic notation and observables of interest, in section 
IIIII we present our simulations and conclusions are given 
in the final section HVl 



II. NOTATION AND PRELIMINARIES 

We simulate the q = 3 Potts model with the energy 
function 
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(1) 



where the sum is over the nearest neighbors of a 3D cubic 
lattice of size L 3 . The spins qi of the system take the 
values qi = 0, . . . , q — 1. The factor of two and the term 
1/q are introduced to match for q = 2 with Ising model 
conventions Q. 

Our simulations are carried out in a multicanonical en- 
semble [8j, covering an inverse temperature range from 
(3 = Anin = (infinite temperature) to (3 = (3 max = 0.33 



below the phase transition temperature. Instead of rely- 
ing on a recursion (see, e.g., (9(), the multicanonical pa- 
rameters were extracted by finite size (FS) extrapolations 
from smaller to larger system, which is efficient when the 
FS behavior is controllable. Frequent excursions into the 
disordered phase all the way to j3 = secure equilibra- 
tion of the configurations around the transition and in 
the ordered phase. 

By the fluctuation-dissipation theorem the specific 
heat is 



C v = «£ 2 ) - (Ef) 



(2) 



We put (3 in parenthesis, because the number 2 is later 
used as a superscript of (3. The energy density per site is 
defined by 



so that the latent heat is 
Ae = 



(3) 



(4) 



where e + and e~ are the energy endpoints of the high (+) 
and low (— ) temperature phase at the transition temper- 
ature T c — 1/(3 C . The entropy density is 



s = (3(e-f) , 



(5) 



where / is the free energy density. As the free energy is 
continuous at the phase transition the entropy gap across 
the phase transition is 



As = (3 C Ae . 



(6) 



To determine the entropy and energy gaps from data on 
finite lattices we follow Ref. [ll| and study the scaling 
of the specific heat maxima C max (L). The leading order 
coefficient 02 of the fit 



C m ax(i) = ai + a 2 L 3 

is related to the latent heat by 

(A) 2 (Ae) 2 
&2 — : 



(7) 



(8) 
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and using ([6]) gives 

As = 2^ . (9) 

The entropy densities in the disordered (+) and ordered 
(— ) phase at the transition temperature are defined by 

s + = Pc (e + - /(A)) , (10) 
a" = A (e- - /(/3 C )) . (11) 

Calculation of the endpoints from Monte Carlo (MC) 
data faces some technical difficulties, which are overcome 
in section UlIB I 

In a MC simulation of the Gibbs canonical ensemble 
the entropy of the system is only determined up to an ad- 
ditive constant, whereas in our multicanonical simulation 
this constant is determined by the known normalization 
at p = 0: 

So = In (3 l3 ) and s = = In 3 ~ 1.098612 . (12) 
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FIG. 1: Histograms of the energy density reweighted to equal 
heights on various lattices. 



III. SIMULATION RESULTS 

In Tabic |T] we list the lattice sizes used, number of 
sweeps (sequential updates of the lattice for which each 
spin is touched once) performed with the multicanonical 
Metropolis algorithm and the number of cycles 

(Pe < /3min) -> (Pe > /3 maX ) -> (Pe < Anin) , 

the Markov process performed during the production 
run, where p e is the effective energy- dependent p of the 
multicanonical procedure. Data are energy histograms, 
which are collected over the statistics given by the second 
number in the production statistics column (sweeps) of 
Table HI This is repeated for 32 bins, the first number in 
this column. The data analysis relies on converting the 
bins into jackknife bins, which are then reweighted to 
canonical ensembles using the logarithmic coding of @. 
Error bars are obtained by performing analysis calcula- 
tions for each jackknife bin. 

TABLE I: Statistics of multicanonical simulations of the 3- 
state Potts model on L 3 lattices. 



L 


sweeps 


cycles 


20 


32 x 1.2 • 10 5 


59 


30 


32 x 5.2 ■ 10 5 


71 


40 


32 x 1.5 • 10 6 


73 


50 


32 x 6.0 ■ 10 6 


131 



A. Interface tension and specific heat maxima 

Strictly speaking there are no phase transitions on fi- 
nite lattices. But one can define pseudo-transition tem- 
peratures Tpt(L) = l/P pt (L), which agree up to 1/L 3 



TABLE II: First inverse pseudo-transition temperature and 
the interface tension on L 3 lattices. 



L 




2cr od 


Cmax 


20 


0.275311 (33) 


0.00152 (14) 


26.27 (51) 


30 


0.275284 (11) 


0.001440 (49) 


62.49 (88) 


40 


0.2752853 (59) 


0.001522 (37) 


136.4 (1.5) 


50 


0.2752838 (23) 


0.001574 (22) 


263.2 (1.7) 



corrections with T c . We give a first definition, pp t (L), as 
the P values which reweight the double-peaked histogram 
at the transition to equal heights. These histograms are 
presented in Fig. [T] and corresponding pz t (L) values are 
given in Table HT1 Another definition, Pp t (L), will be in- 
troduced later. For the first-order phase transition the 
FS behavior of these pseudo-transition definitions is (up 
to higher orders in l/L 3 ) 

p; t (L) = pi + ^, i = l,2. (13) 

Fitting the Pp t (L) values of Table |TT] to this equation 
yields 

Pi = 0.2752827 (29), b 1 = 0.15 (22) (14) 

with a goodness of fit Q Q — 0.89. Within statistical 
errors the coefficient b 1 is zero, indicating that the FS 
corrections to P\ are mild. 

For L —> oo the interface tension between ordered and 
disordered phases is [l(| 

where P max (L) represents the value of the maxima when 
the energy histogram is reweighted to equal heights and 
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FIG. 2: Entropy density of the 3D 3-state Potts model on 
different lattices. Horizontal lines represent s + and s~ values. 
FS effects are not visible on the scale of this figure. 
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FIG. 3: Entropy density in the vicinity of the phase transition. 
The height of the bar on the right represents the entropy gap 
As. The horizontal lines are s + and s~ and the vertical line 
is/3 c 2 . 



Pmin(L) the minimum in between the peaks (Fig.Q}. The 
results are contained in Table [TTJ A fit to the form 

2a od {L) = 2a od + -^ (16) 

gives 

2a od = 0.001602 (35) , c= -0.100 (52) (17) 

with Q = 0.22. This is consistent with previous lit- 
erature: 2a od = 0.001568(52) of @ (Q = 0.59) and 
2<j od = 0.00163 (2) of @ (0 = 0.49), where the Q values 
are now from Gaussian difference tests |9j with our fit. 
However, if one includes capillary waves [12l. [l3l. [l4| the 
fit form becomes, see Eq. (16) of [Tsl ]. 

2a od (L)+ ] ^-=2a od + ^ (18) 

and we obtain a different consistent estimate: 

2a od = 0.001806 (35) , c= 1.379 (52) (19) 

with Q = 0.96. 

Table |H] collects also our estimates of the maxima of 
the specific heat C max (L). Fitting them all to Eq. ([7|) 
gives the rather low goodness of fit Q = 0.028. Following 
standard practice we omit therefore the smallest (L = 20) 
lattice from the fit and find 

ai = 7.0(1.3) and a 2 = 0.002044 (20) (20) 

with Q = 0.28. 

B. Critical temperature, entropy and energy 

The entropy density from our multicanonical simula- 
tions is shown in Fig. [5] Error bars and FS corrections 



are not visible on the scale of this figure and the vicinity 
of the phase transition is enlarged in Fig. [3] 

To determine the endpoints of the entropy and energy 
on the disordered and ordered side of the transition, we 
introduce a second definition of inverse pseudo-transition 
temperatures, j3p t (L), which differ from the equal heights 
definition in crucial details. We adjust the weights of the 
double peak histograms so that 

eL(P 2 pt ) = l[et(f3 2 pt )+e- L ((3 2 pt )} (21) 

holds, where the central energy density ex,(/3^ t ) is the ex- 
pectation value at 0p t (L) and e L (f3p\) are the locations 
of the maxima of the double peak histogram at 0i t (L), 
e\ on the disordered and on the ordered side of the 
transition. The idea behind this construction is to ensure 
that the energy endpoints ^ are positioned symmetri- 
cally about the central energy density: 

ef=e L (p 2 c )±±Ae L . (22) 

In the following we call this procedure equal distances 
reweighting. Using ([5]) and ^ one finds that 

sf = s L (f3 2 c )±lA SL (23) 

holds as well. 

The histograms reweighted to equal distances are 
shown in Fig. 2] and the corresponding f3p t (L) are col- 
lected in Table IIII1 For the L — 20 histogram of the fig- 
ure (3 — 0.275311 is used just to illustrate that the double 
peak disappears before condition (f2T|) gets fulfilled. This 
is the reason why there are no L = 20 values in the table. 
Fitting the (3 2 C (L) values of Table HH to (13]) yields 

f3 2 = 0.2752724 (44), b 2 = -2.61 (40) (24) 
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FIG. 4: Histograms of the energy density reweighted to equal 
distances on various lattices. 



TABLE III: Second inverse pseudo-transition temperature, 
central energy and and entropy densities on L 3 lattices. 



L 








h{P 2 pt ) 


30 


0.275166 (13) 


-1.336 (14) 


0.8522 (37) 


-4.43275 (14) 


40 


0.2752369 (52) 


-1.345 (11) 


0.8497 (29) 


-4.431840 (62) 


50 


0.2752508 (22) 


-1.3444 (51) 


0.8498 (14) 


-4.431652 (26) 



with a goodness of fit Q = 0.21. The estimates (|14j) and 
((24)) are at the edge of being consistent (Q = 0.051 for 
the Gaussian difference test). Our second estimate for 
(3 C is less accurate than our first, apparently because the 
location of the newly defined maxima of the energy prob- 
ability density is less stable than for the equal heights def- 
inition. Also FS corrections to the second estimate are 
not negligible as those to the first. In absolute numbers 
the difference is less than 1 in the fifth significant digits, 
so it does not propagate into the (3 C factor used for the 
far less accurate latent heat. Still the second definition 
is needed to enable estimates of e ± and s . 

To give a combined estimate of j3 c , we average the 13], 
and flc estimates weighted by their inverse variances and 
keep the error bar [16| of the more accurate fi\ result to 
obtain 



(3 C = 0.2752796 (29) . 



(25) 



For comparison, the weighted average of on lattices 
L = 10, . . . , 36 is /3 C = 0.2752825(50) and'the Gaussian 
difference test with our value gives Q = 0.62. 

To calculate As and Ae we define jackknife estimators 
on finite lattices for them: 



As L 



C max (.Zy) 



L 3 



(26) 



Ae, 



L 3 



(27) 



Their values will be needed for combination with jack- 
knifed central entropy and energy values, so error bars 
for and s ± can be obtained. Only L > 30 lattices 
are used as we have already seen for that the smallest 
lattice spoils the fit. Estimates of the means and their 
error bars follow in the usual jackknife way. Fitting them 
with 1/L 3 corrections, the entropy gap is 



As = 0.09045 (41) 



(28) 



It is represented by the height of the filled bar on the 
right of Fig. [3l In the same way we find for the latent 
heat 



As 



= 0.3286(15) 



(29) 



This value is somewhat larger than Ae = 0.32320 (94) of 
(Q = 0.002) and in agreement with Ae = 0.3291 (33) 
of 0] where additional data of [3| on lattices up to L = 
48 was taken into account (Q = 0.89). However, these 
results depend still beyond the range of one error bar on 
the order in which the functions are evaluated. Using (J9j) 
and ijHJl to estimate the gaps from a-i of Eq. (f2"0"|) , one finds 
As = 0.08933 (38) and Ae = 0.3245 (14). This Ae value 
is consistent with both results of the previous literature 
as it is in the middle of them. 

The central energy (|21[) . entropy and free energy den- 
sity values (all defined at 0L) are listed in Table IIIII 
From 



e L (P 2 pt ) 



(30) 



we find the infinite volume extrapolation 

e{0t) = -1.3470(74), Q = 0.82. (31) 

Similarly we determine 

s{P 2 c ) = 0.8491 (21) , Q = 0.82, (32) 

and the corresponding (J5]) free energy density 

f(Pl) = -4.431364 (50) , Q = 0.18. (33) 

The statistical errors of ftiPpt) are quite small due to 
correlations between ei(/3p t ) and Si(/3p t ). 

From (12"2")) and (f2"3"|) we find finite lattice estimators ej 
and s^ , which are collected in Table HVl Fits with 1 / L 3 
corrections give: 



e+ 


= -1.1826(73), 


Q = 0.74, 


(34) 


e~ 


= -1.5112(79), 


Q = 0.88, 


(35) 


s + 


= 0.8943 (21), 


Q = 0.73, 


(36) 


s~ 


= 0.8038 (22), 


Q = 0.90. 


(37) 



The values of s + and s~ are shown in Figs. and [3] as 
horizontal lines. The vertical line in Fig. [3] represents the 
inverse transition temperature 
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TABLE IV: Results for J22J and (23]). 



L 


4 






S L 


30 


-1.161 (13) 


-1.511 (15) 


0.9003 (36) 


0.8041 (39) 


40 


-1.177 (10) 


-1.513 (12) 


0.8958 (28) 


0.8035 (31) 


50 


-1.1777 (50) 


-1.5111 (53) 


0.8957 (14) 


0.8039 (15) 



IV. SUMMARY AND CONCLUSIONS 

The main result of this work is the entropy density with 
proper normalization, shown in Fig. [2 where sq = In 3. 
The transition region is enlarged in Fig. [3] The values 
of the entropy density s + (|36p on the disordered and s~ 
(13"?)) on the ordered side of the transition in % of sq are, 



respectively, 82% and 73%. Thus, there are 3 states per 
site at infinite temperature and effectively 2.45 states on 
the high- and 2.24 states on the low-temperature side of 
the phase transition. Our corresponding estimates are 
given in Eqs. (J34J) and ([35| . 

Other results are (3 C (|25|) , the free energy density at the 
critical point (|33|) . the entropy gap (|28|) . the latent heat 
(|2U)) and the interface tension (fT9"| . It is notable that the 
inclusion of capillary waves enhances the estimate of the 
interface tension by more than 10%. 
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